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Abstract 

We consider the system 

— Auj + a(x)uj = HjV?j + (3 ^ u\uj, uj > 0, j = 1,... , n, 



on a possibly unbounded domain Q C M. , N < 3, with Dirichlet boundary 
conditions. The system appears in nonlinear optics and in the analysis of 
mixtures of Bose-Einstein condensates. Two components Ui,uj are said to 
be locked if Ui/uj is constant. The main results are concerned with the bi- 
furcation of solutions (ui, . . . , u n ) where some components are locked from a 
branch where all solutions are locked. The bifurcation parameter is (3. 
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1 Introduction 

The system of coupled Gross-Pitaevskii equations 

(1.1) - id t 7pj(x, t) = Axtpj + njltpj^ipj + (3 ^ iV'fcl 2 ^ 3 = 1, • • • , n, 

in M N , N < 3, with parameters fj,i, . . . , fi n > 0, (3 e K, has found considerable 
interest in the last years. It appears in nonlinear optics and in models for mixtures of 
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Bose-Einstein condensates; see [1, 7, 15] for physics, and [2, 4, 5, 6, 10, 11, 14, 17, 18] 
for mathematics papers. We deal with the case of attractive self-interaction, i. e. 
Hi > 0, and repelling interaction between different components, i. e. (3 < 0. The 
ansatz 

il>j(x,t) = e %t Uj(x), j = l,...,n, 

for stationary waves leads to the following elliptic system for the amplitudes 
Mi, • • • ,u n : 



We keep /ii, . . . , \x n > fixed and take (3 as bifurcation parameter. 

Most of the above mentioned papers papers deal with the case n = 2 of two 
components. Using variational methods, the existence of ground and bound states 
is obtained. In [4] a different approach is pursued using bifurcation methods, but 
only for n — 2. We shall extend this bifurcation approach from n = 2 to arbitrary 
n > 2. This leads to interesting new features and difficulties. For instance, in the 
radial setting of [4] the linearizations have a one-dimensional kernel at the bifurcation 
parameter. This allows to use degree theory to prove global bifurcation of solutions. 
For more than two equations we have to deal with high-dimensional kernels forced by 
the structure of the system. We shall partially explain this with a hidden symmetry 
of the problem which will lead to multiple global bifurcation branches. 

Actually, we shall deal with a more general problem. Let Q C M. N , N < 3, 
be a domain which is invariant under a closed subgroup G C O(N), and suppose 
a G L^ C (Q) is invariant under G. We require that —A + a is a positive self-adjoint 
operator on L 2 (Q), and define 



We also require that the embedding E •— > L A (Q) is compact. This is the case, 
for instance, if Q is bounded and G is trivial, or if Q is radially symmetric and 
G = 0(N), 2 < N < 3. It is also the case if Q = R N and a(x) ->• oo as \x\ ->■ oo. 



(1.2) 




£:={u6 : u is G-invariant, J Q au 2 < oo} . 



The norm in E is given by 
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The problem we investigate is 

{- Auj + a(x)uj = fiju* + (3 u l u i » j = 1, . . . , n, 

G -E, Mj > 0. 

As mentioned above, this covers the important special case Q. = M. N , N = 2 or 
N = 3,E = H^ ad (R N ). 

A solution (ui, . . . ,u n ) of (1.3) is said to be locked if Ui/uj is constant for all 
In this paper we first describe the set of locked solutions of (1.3). We then 
investigate the bifurcation of non-locked solutions from the set of locked solutions. 
In particular, we shall prove the bifurcation of partially locked solutions where Ui/uj 
is constant for some i ^ j, but not for all indices. 

Without loss of generality we assume < jii < /x 2 < • • • < fi n throughout the 
paper. 



2 Branches of locked solutions 

Given a solution of the scalar equation 
(2.1) 



— Aw + a(x)w = w 3 
w e E, w > 



there exists a branch T w C M x E n of locked solutions (/3,ui, . . . ,u n ) of (1.3) such 
that each Ui is a multiple of w. In order to describe this branch we set 

(2.2) g{fi ).. = i + fij2-}— 
and for j = 1, . . . , n: 

(2.3) aj (J3) := {fa - 0)g(J3))- 1/2 

Observe that g is defined and strictly increasing in the interval (— oo,/ii), and that 
it satisfies g(0) = 1, g(j3) — > 1 — n < as (3 — > — oo. Thus there exists a unique 
(3 = j3(fj,i, . . . , fj, n ) < such that g(f3) = 0. Consequently, the functions ai, . . . , ot n 
are defined for j3 G (j3, fXi) U (/i n , oo). 

Proposition 2.1. Problem (1.3) has a locked solution (/3, Ui, . . . , u n ) precisely in 
the following cases: 
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(i) (3 G /xi) U (fi n , oo): Then Uj = aj((3)w,j — 1, . . . ,n, for some solution w of 
(2-1). 

(ii) (3 = fii = . . . = fi n : Then Uj = otjW, j — 1, . . . , n for some solution w of (2.1), 
ai,..., a n > 0, a\ + . . . + a n = 1/(5. 

Proof. If (/3, ui, . . . ,u n ) is a locked solution then u\ solves — Ami + a(x)u\ = cu\ 
for some constant c > 0. It follows that w = c 1 / 2 -^ solves (2.1), and all u/s are 
multiples of w. 

Now write Uj = ajw, j = 1, . . . ,n, where w is a solution of (2.1). Then (1.3) 
leads to 



cnjiv 3 = (—A + a(x))ajW 



hence, 



(2.4) Hjtf + p^a^l, j = l,...,n. 
This implies 

(2.5) (fii - (3)a* = (fij - f3)a] for all i, j. 

1 /2 

In case (i) it follows from (2.5) that a,- = [j^z^j «i for j = 1, . . . , n. A simple 
computation using (2.4) now yields 

ax = (0/i- f3)g(/3))- 1/2 = ai ((3), 

and (Xj = otj((3) for all j follows immediately. This proves in case (i) that a locked 
solution (j3, Mi, ... , u n ) necessarily has the form Uj = a.j(j3)w for some solution w of 
(2.1). If (i) does not apply then we must have j3 = «i = . . . = « n , which implies 
f3(al + . . . + a n ) = 1, again by (2.4). This is as claimed in (ii). 

Finally, an elementary calculation shows that (f3,U\, . . . ,u n ) as described in (i) 
or (ii) is a solution of (1.3). □ 

Now we fix a solution w oi (2.1) and set 

u(l3) := (a 1 (f3)w, a n ((3)w) for f3 G (f3, «i) U («„, oo). 

We want to investigate the bifurcation of solutions of (1.3) from 

T w := {((3, u((3)): (3 £0,^)13(^,00)}. 
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3 Necessary conditions for bifurcation 



We need to linearize (1.3) at (P,u(P)) G T w . A simple calculation leads to the 
system 



(3.1) 
with 



(-A + a(x))<f> = w 2 C(/3)(f) 
= (0 1 ,...,0„) T GE" 



C(P) = E n + -^D((3) G R nxn . 



Here E n is the n x n unit matrix and 



( Hill Pl\l2 ••• Pllln\ 
Mll2 Vlll ••• Pl2ln 



\hlln Pl2ln ■■■ Hull/ 

where ^ = -fj(p) := (/ij - f3)~ 1/2 . Hence, D(p) = (dij(P))ij=i,..., n with 

P 



({jH-P^j-P)) 



1/2 



Hi 



Hi- p 

Observe that C(p) and D(P) are symmetric matrices. 
Lemma 3.1. a) C(P) has the eigenvalues 3 and 

2 2 



if i ^ j\ 
if % = j. 



f(P) ■= 1 + 



i + 



b) If P ^ £/ien i/ie eigenvalue 3 simple with eigenvector 

(3.2) 6!(/3) = ( 7 i 09), ■ ■ ■ , 7n(/?)) T = ((^i - /T 1/2 , ...,(*>- P)- l,2 )\ 

and the eigenvalue f(P) has multiplicity n — 1 with eigenspace bi(P) ± . If P = 
then /(0) = 3 and C(0) = 3£„. 
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Proof. A simple calculation shows that D(f3)bi(f3) = g(/3)b 1 (P), hence C(f3)bi((3) = 
3bi(P). Moreover, setting 

r 7 ,(/3) ak = i, 

(3.3) bj((3) = (bj!, b jn y where b jk = < -71 (/3) if k = j, 

[0 else. 

for j = 2,...,n we have 6j(/3) JL h (f3) and D(p)bj(p) = bj(p). This implies 
C(/3)6 J -(/3) = /(/9)6 j (/3)forj = 2,...,n. 

Since /(/?) 7^ 3 for j3 7^ we see that C(j3) has the eigenvalue f(/3) with 
eigenspace span{6 2 , • • • , b n } = 6i(/3)- L . If j3 — then /(0) = 3, hence 3 is the 
only eigenvalue of the symmetric matrix C(0). □ 

Remark 3.2. Observe that f : (/3,/ii) — >■ (1, 00) is a strictly decreasing diffeomor- 
phism. In fact, 

Moreover, f(/3) ->■ 1 as /3 /ii, and /(/3) ->■ 00 as /3 \ (3. 

Remark 3.3. Since C{j3) is a symmetric matrix depending smoothly on (3, and 
since the eigenvector bi(j3) also depends smoothly on (3, there exists a smooth map 
T : (/3,/ii) ->• SO(n) such that 

T(P)- l C{P)T(P) = diag(3, f(J3), /(/?)). 

A point (f3,u(f3)) e T w can be a bifurcation point of solutions of (1.3) only if 
(3.1) has a nontrivial solution 0. Nontrivial solutions of (3.1) are closely related to 
the weighted eigenvalue problem 

(3.4) - Aip + a(x)^ = \w 2 tp, ip G E. 

Remark 3.4. Recall thatw is a non- degenerate solution of (2.1) if and only if X = 3 
is not an eigenvalue of (3.4). In the degenerate case where ip 7^ is a solution of 
(3.4) with X = 3, we see that <f) = b x {(3)ip G E n , bi((3) G W 1 as in (3.2), is a solution 
(3.1) for every j3 < j3 < fj,i. If w is not an isolated solution of (2.1), say uik — >■ w 
for a sequence of solutions of (2.1), then we have branches T Wk of solutions of (1.3). 
In this case every point on T w is a bifurcation point for (1.3). The problem is more 
subtle if w is an isolated, degenerate solution of (2.1). This case will not be treated 
here. 
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From now on we assume that w is a non-degenerate solution of (2.1). Then (3.4) 
has a sequence of eigenvalues Ai = 1 < A 2 < A 3 < . . . such that \ k ^ 3 for all k G N 
and Afe — > oo as k — > oo. The multiplicity of A& as eigenvalue of (3.4) is denoted by 

(3.5) nfc = dim\4 with Vk = kern(— A + a — XkW 2 ). 

Proposition 3.5. a) A point ((3,u((3)) e T w is a bifurcation point only if (3 = 
(3k := f~ l (\k) for some k > 2 with f from Lemma 3.1. 

b) If (3 = (3k then <f> = (0i,...,0„) G E n solves (3.1) if, and only if, <pj G 
Vk for j = l,...,n, and _L b\{(3) with b\{(3) G R™ from (3.2), that is, 
Sj=i 7j(/3)<Aj = 0. In particular, the space of solutions of (3.1) has dimension 
(n - 1) • n k . 

Proof, a) Let T((3) G S'O(n) be as in Remark 3.3. Then solves (3.1) if, and only 
if, ip = T(/3)" 1 solves 

f -A^i + a^i = 3w 2 fa 
1 ' J \ -A^- + afa = f((3)w 2 fa, j = 2, . . . , n. 

Since w is non-degenerate this implies fa — 0, hence a nontrivial solution -0 exists, 
if and only if f((3) = \ for some k > 2. (The case fc = 1, i. e. Ai = 1, (3 = fxi does 
not apply; see also Remark 3.6.) If f((3) = \ k then ip = (fa, . . . ,fa) T with *0i = 0, 
V'fc, solves (3.6). 

b) We fix (3 = f~ 1 (Xk) and write T = T((3) = (T 1? . . . , T„). The transforma- 
tions ip i — ^ T-0, h-> T^ 1 ^, are isomorphisms between the solutions of (3.6) and 
(3.1). As a consequence of Lemma 3.1 we see that T\ is a multiple of &i(/3), and 
span{T 2 , . . . ,T„} = &i(/3) x . Thus solves (3.1) if, and only if, ip = T~ l <p = T T (f) 
solves 

(3.7) 0) = fa = 0, hence = (6i(/3), 0) = 7 i Wi + • • • + 7«(/3)0„ 
and 

(3.8) (T J ,<j ) )=faeVk, for j = 2,...,n. 

Recall that the eigenvectors 6 2 , ■ ■ ■ ,b n from (3.3) satisfy span{6 2 , • • • , b n } = 
span{T 2 , . . . ,T n }. It follows that (3.8) is equivalent to 

(3.9) (&;,0>=7iO0)0i-7i(/%e^ for j = 2,...,n. 
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Now (3.7) and (3.9) are equivalent to 



n 

0i, . . . , <p n G V k and ^ lj{P)4>j = 0. 

i=i 

□ 

Remark 3.6. By our definition T w does not contain a point with the parameter value 
f3 = n\. If Hi < fi n then it is not difficult to show that (/i l5 u±, . . . , u n ) G T w implies 
that at least one Uj = 0. // fii = fi n then (yU 1? (n/ii) _1 / 2 u>, . . . , (n/ii) _1 / 2 u>) G T w . 
This is a bifurcation point where the bifurcating locked solutions are described in 
Proposition 2.1: 

n 

(fj,i, aiw, . . . , a n w) with aj > 0, a 2 = 1/(3 = \ j H\. 

3=1 

4 Sufficient conditions for bifurcation 

We first recall the variational structure of problem (1.3). Since N <3 the functional 
Jp(ui, ...,u n ) 

=\£[ a wi 2 + «(*)«?) -\±f -fx;/ 

3=1 jQ 3=1 Jn i<j jQ 

1 r l "\ r r 

= 2^\WA\l-i^»j\\ u j\\i*-2l2 / u ^ 

3=1 3=1 i<j 

is well defined for u\, . . . , u n G E. It is well known that Jg : F" — > R is of class C 2 , 
and critical points of Jg are (weak) solutions of (1.3). 

For (3 G (/3,/ii) U (/x n ,oo) let m(/3) be the Morse index of w(/3) = 
(aii (/?)«;, . . . , a n (j3)w) as critical point of J/?. Recall the bifurcation parameters 
from Proposition 3.5. 

Lemma 4.1. For k > 2 i/iere /ioWs 

m(/3 fe — e) — m(/3 fc + e) = (n — l)n fe /or e > small. 
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Proof. Consider the quadratic form on E n given by 

■ =1 Jn Jn 

= u«»- [ w\cm^) 

Jn 

where C(j3) is as in (3.1). Then m((3) is precisely the index of Qp. Let E n = 
Vp k ® V® k © Vg~ be the decomposition of E n into the positive eigenspace Vp of Q/? fe , 
the negative eigenspace V g ~ of Qp k , and the kernel 

(4.1) Vp° k = l^eE n : 0, G \4 for j = 1, . . . ,n, = j . 

The derivative Q'p = -§pQp of Qp with respect to f3 is given by 

W) = - / ^ 2 (cw,0)- 
in 

Recall that C(/3) is a smooth function of /3. Now we have 

Qp = Qp k + 03 - + O (I0 " &l) as £ -> 

This implies that Qp > on Vp and Q/j < on Vg~ , if /3 is close to (3^. Thus the 
lemma follows if Q'p k is positive on Vp k . In order to prove this let T(j3) G SO(n) be 
as in 3.3, that is, T depends smoothly on (3, and satisfies 

T(/3)- 1 C(/3)T(/3) = diag(3, f(J3), . . . , /(/?)) =: C T {fi). 

It follows that C(f3)T(f3) = T(/3)C T (/3), hence 

C"(/3) = T\P)C T (P)T{P)- 1 + T(/3)^(/3)T(/3)- 1 - C(/3)T'( j 9)r( j 9)- 1 

For an eigenvector u G &i(/3) _L C lR n of C(/3) corresponding to the eigenvalue f(/3) 
we compute 

T'((3)C T ((3)T((3)- l u = f({3)T'(P)T(P)- l u, 



T^CMT^u = f'(P)u, 
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and 

(C{p)T'(p)T{p)- 1 u,u) = {T\(5)T{(5)- l u,C{(5)u) 

= f((3){T>((3)T((3)-\u). 

This implies 

(C'(f3)u,u)=f(/3)\u\ 2 . 
According to Proposition 3.5 a function G Vg satisfies 

n 

^7j(/3fc)0j = 0, i. e. 0(x) G ^(/^fc)" 1 C IR n for every x eVt. 

3=1 

Here we used Proposition 3.5 b). Consequently, for G Vg we have 

= - / ™ 2 <C"(/3 fe )0,0> = -/'(&) / |0| V > 

Jn Jn 

because f'(Pk) < by Remark 3.2. Therefore Q'^ is positive definite on Vg fc , as 
required. □ 

Theorem 4.2. For every k > 2, the point (j3,u(/3)) G T w is a bifurcation point of 
solutions of (1.3). 

Proof. We can apply a classical bifurcation theorem of Krasnoselski [9] in the version 
of [16, Theorem 8.9]. The argument proceeds as in [4, pp. 354-355]. Alternatively 
one can use [3]. Both arguments yield that the bifurcating solutions are actually 
positive. □ 

Remark 4.3. a) If (n — l)n^ is odd then there exists a connected set C M.xE n \T w 
of solutions (P,u) of (1.3) such that (/3k,u(/3k)) G Sk- This can be proved using 
degree theory in the spirit of Rabinowitz' global bifurcation theorem. Whether Sk is 
unbounded in the (3-direction, or in the E n -direction, or whether it returns to T w is 
unclear in general. 

b) The case n = 2 has been treated in [4j. Clearly, if n is odd degree theoretic 
methods do not apply. On the other hand, one may hope for an (n-l)nk-dimensional 
manifold of bifurcating solutions, or multiple bifurcating branches. In many cases, 
a high dimensional kernel is not generic but forced by some symmetry. In the next 
section we present some hidden symmetry, which is not inherent to the full func- 
tional, but to a special type of solutions. In this way we obtain multiple bifurcating 
branches. 
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5 Branches of partially locked solutions 

In this section we investigate partially locked solutions of (1.3), i. e. solutions (/3,u) 
where Ui/uj is constant for some, but not all i ^ j. If Ui/uj is constant we say that 
Ui and Uj are locked. 

Lemma 5.1. If (j3,u) solves (1.3) andui,Uj are locked then 



Uj = Ui = — — Ui . 

Proof. We set Uj = •yui for some 7 > 0. Then the equations (1.3) for Ui and Uj lead 
to 

-Aui + aui = (fii + /?7 2 )m- +13^2 u ^ Ui 

and, respectively, to 

n 

—Aui + aui = (/ij7 2 + (5)u\ + (5 ^ u\ui 

It follows that 

>-£Y /2 im 
7 - . - 



□ 

The next result reveals some hidden symmetry in the problem which will allow 
a dimension reduction argument below. It is simple to prove but has not been 
observed before. 

Lemma 5.2. Suppose (j3,u) solves (1.3) for some i G {1, . . . ,n}, and suppose 

'^-PV'\ _7,(/3) 



Uj = I I Ui = Ui 



Then (1.3) also holds for Uj. 



Proof. We set 7 := ^Hy, so Uj = 7^. Multiplying the equation for Ui 



-A^ + aui = Hi/if + 13^2 u l u i = + I3l 2 )u 3 i + /3 

k+i k^i,j 



U k Ui 
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by 7 yields 

-A Uj + auj = + M 2 )iu\ + u k u i = (Vjl 2 + Phuf + P Yl u ^ u i 

= n jU ) + {3j2u 2 kUj . 

This is equation (1.3) for Uj as claimed. □ 

Let V = {Pi, . . . ,P m } be a partition of {1, . . . , n}, i. e. {1, . . . , n} — PiU. . .UP m 
is a disjoint union and Pi, ... , P rn ^ 0. We write \V\ := m for the cardinality of V . 
For (3 < we define 

Xj := G E n | \/k = 1, . . . ,m V i, j G P fe : Uj = ■ 

A solution -u G Xj of (1.3) will be called P- locked. 

Proposition 5.3. If u G Xj is a critical point of Jg|Xj £/ien u is a critical point 
of J/3, hence a solution of (1.3). 

Proof We may assume ttj G p for i = 1, . . . , m. That w G Xj is a critical point of 
Jp\X~p is equivalent to tii, . . . , u m being solutions of (1.3) with Uj = ^j^Ui j G p, 
i = 1, ... ,m. Lemma 5.2 implies that w = (wi, . . . ,«„) solves (1.3), hence u is a 
critical point of Jg. □ 

One can now use Proposition 5.3 to find bifurcations of partially locked solutions 
in Xj. We may assume that i G Pi for i = 1, . . . , m. Define 



i 



■ P ■ E m ->■ Xj, (Ui, . . . , Um) H> ...,«„), 



by 



M i = ^4^ u i if J e P, i = 1, . . . , m. 



Let Ti : P n — > P m be the projection onto the first m components, so that 7r o is 
the identity on E m for every j3 < fii. Finally we define 



: = J o i? : P m -)• 
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for (5 G R. A critical point u G E m of yields a critical point ip(u) G Xj of 
Jp\X£, hence a P-locked solution of (1.3) as a consequence of Proposition 5.3. The 
branch 

-C = {(/3,nu) : (/3,n) G T w } = : (&#(«)) G T w } 

is a branch of critical points corresponding to the branch T w of locked solutions. 

Let m v ((3) be the Morse index of nu(/3) = {ol\{0)w, . . . , a m (/3)w) as critical point 
of Jj. Bifurcation of P-locked solutions occurs if m v (j3) changes. 

Lemma 5.4. For k > 2 t/iere ZioZds 

m V (/3k — e ) — mV (l3k + e) = (|P| — / or £ > small. 

Proof. This follows from Proposition 5.3 and Lemma 4.1 if we can show that 

dim(V£ nX£) = (|P|-l)n fc ; 

here Vg is the kernel of the Hessian of Jp k at n(/3 fc ) as in (4.1), and = dim 14 
as in (3.5); recall the description of V® k in Proposition 3.5. For simplicity we write 
lj = lj(Pk) below. 

n 

</>eVjj k nx% k <=^> 0i, . . . , <p n g v k , = °> and 

0o = — & for j £ Pi, i = 1, . . . , to = P 

7?' 



)!,..., m G Vfc, 53 I — J = 0, aIld 



= ij fc (^l, ■■■Am) 



□ 



Corollary 5.5. For every partition V of {1, ... ,n} with |P| > 2, and for every 
k > 2 i/ie point E T w is a bifurcation point of V -locked solutions of (1.3). 

It is a global bifurcation point of V -locked solutions if (|P| — 1) ■ n& is odd. 

In the application of the Corollary 5.5 one has to be cautious because a P-locked 
solution may also be P'-locked for P ^ V. This is the case if for any P[ G V there 
exists a Pj G V with P[ C P,. On the other hand, the bifurcating P-locked solutions 
are different from the P'-locked solutions if every solution which is both P-locked 
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and P'-locked is automatically completely locked. Recall that the completely locked 
solutions are those in T w . 

For a subset Ad {!,..., n} we set A c := {1, . . . ,n} \ A and V A ■= {A, A c }. 
Given two subsets ^ A, B C {1, . . . , n) such that A ^ B and A c ^ B, a solution 
which is "P^-locked and "P^-locked is necessarily completely locked. We therefore 
obtain the following multiplicity result. 

Corollary 5.6. Suppose is odd. For every subset ^ A C {1, . . . , n} there exists 
a global branch oj "V \-locked solutions of (1.3) bifurcating from T w at u(j3k))- 
The branches and are disjoint except when A = B or A = B c . In particular, 
there exist at least 2 n ~ 1 — 1 such global branches which are different. 

Remark 5.7. a) The corollary applies, for instance, if'QcWisa bounded interval, 
or ifVL C 1^, 2 < N < 3, is radially symmetric and one looks for radially symmetric 
solutions, i. e. G = O(N). Then problem (2.1) has a unique solution which is also 
nondegenerate in E. This has been proved in [12, 13, 19, 8], in the case of a ball, 
an annulus, and M. N . 

b) The global features of the branches are not known. In the radial setting 
and n = 2 it has been proved in [4] that the bifurcating branches Sk = C 
(—oo,0) x E 2 are bounded over bounded subsets of {— oo,0) ; hence they cover all of 
(— oo, f3k). An extension of this result to n > 2 is work in progress. 

Acknowledgment: The author thanks Pavol Quittner for discussions on the 
topic, and for the invitation and hospitality during a visit at the Comenius University 
in Bratislava in spring 2011. 
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